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SYMBOLS 


a acceleration of flow 

A variable used to define equation (3.29) 

B variable used to define equation (3.23) 

2 

c speed of sound, c = yp/p 

Cg speed of sound at x = 0, t = 0 

C” characteristic line with dx/dt = u - c 

C + characteristic line with dx/dt = u + c 

D variable used to define equation (2.32) 

E variable used to define equation (2.33) 

h enthalpy 

k(n) acceleration of a characteristic line 

m mass flow rate 

M Mach number 

OA head characteristic of the expansion fan 

OB tail characteristic of the expansion fan 

p pressure 

q energy formation 

Q flow properties 

R universal gas constant Rq divided by the effective molecular 

weight of the particular gas 

T temperature 

u velocity of flow 

v relative velocity 

x displacement 

a dimensionless heat of reaction, a = 


iii 



y 


T1 

5 

5 2 

TT 

P 

V 


ratio of specific heat, y = c^/c^ 
angular coordinate, q = x/c^t 

slope of the characteristic line near the origin 

slope of the tail characteristic line OB at origin 

slope of the head characteristic line OA at origin 

dimensionless final pressure ratio, P 1 /P 2 
density 

specific volume ratio, ^2^1 


Subscripts 

d reaction front 

p piston 

s shock 

t partial derivative with respect to time 

x partial derivative with respect to space 

1 region behind the shock wave or region behind the expansion fan 

2 region ahead of the wave, i.e., free stream 

3 region inside the expansion fan 

4 region between the Chapman-Jouguet detonation and expansion fan 

5 region between precompression wave and Chapman-Jouguet deflagration 

wave 


Superscripts 

(0) leading term 

(1) first order term 

- normalized quantity 

~ flow properties on higher pressure side 

* Chapman-Jouguet reaction front 
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CHAPTER 1 


INTRODUCTION 


1 . 1 Purpose of Paper 

This paper presents derivations of canonical solutions for the initial 
value problem of one-dimensional unsteady gas dynamics with discontinuous 
nonuniform initial data. The canonical solutions are constructed by assuming 
(i) regular power series expansions for regions ahead of and behind each 
discontinuity wave, namely, a shock wave (with or without chemical reaction) 
or a contact line, and (li) special expansions within a centered expansion 
fan. The solutions define not only the speeds but also the accelerations of 
the discontinuity waves and the expansion waves. These canonical solutions 
are intended to be employed as the building blocks in a numerical scheme for 
the solution of a general initial value problem. To understand the need for, 
and the importance of, this work, a brief review of related methods is 
included in the next section. 


1.2 Review of Work 

A number of fluid dynamics problems as well as other physical problems 
are dominated by the presence of strong discontinuities. In order to compute 
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the solutions to these problems accurately, the discontinuities have to be 
understood and specific information about their generation and propagation 
must be included in numerical solution schemes. Many computational methods 
for gas flow are based on approximating the problem with a number of 
elementary or canonical flow problems called Riemann problems. An example of 
a Riemann problem is the simple shock tube problem for which initial data, 
discontinuous but otherwise uniform, is given. The canonical problem provides 
an explicit and elementary class of solutions which contain extensive 
information about wave interaction. The existence and uniqueness of the 
solution of the Riemann problem for gas dynamics, subject to an appropriate 
formulation of the entropy condition, was established by Liu [1]. Glimm 
recently developed a front tracking code with the aim of providing a general 
and flexible method for obtaining accurate solutions to problems which are 
piecewise smooth [2], The Riemann problem provides the key input to this 
method. Chorin introduced the random choice method as a numerical tool for 
solving hyperbolic systems [3]. The method is based on a constructive 
existence proof due to Glimm [4]. The solution is advanced in time by a 
sequence of operations which includes the solution of Riemann problems and a 
sampling procedure. Chorin further illustrated the usefulness of the method 
in the analysis of reacting gas flow [5]. Examples were given of time- 
dependent detonation and deflagration waves, with infinite and finite reaction 
rates. 

For a steady supersonic flow, the flow quantities and their gradients 
behind a shock can be related to the flow quantities ahead of the shock and to 
the slope and the curvature of the shock [6]. These relationships can be used 
in the numerical computation of steady flow fields with shock waves, and 
because they explicitly account for nonuniformities on either side of the 



shocks, they are inherently superior to uniform flow relations. A more 
extensive application of these ideas was carried out by Darden [7] in 
conjunction with calculation of sonic boom in the vertical plane of symmetry 
near a point of coalescence of two shocks. Near such a point, there exist 
regions bounded by shocks or by a shock and a slip stream or an expansion 
wave. To include asymmetric lifting effects, the second circumferential 
derivatives of the flow quantities in each region in the vertical plane are 
needed. Explicit relationships between the various derivatives of the flow 
quantities at the point of coalescence and the slopes and curvatures of the 
shocks, the slip streams, and the expansion fans were derived in reference 7. 

In the current paper, canonical solutions for one-dimensional gas 
dynamics problems with discontinuous and nonuniform initial data are 
presented. They are valid in a small space-time region around the 
discontinuities. The canonical solutions relate the first partial derivatives 
of the flow quantities with respect to t and x to the velocities and 
accelerations of the resulting discontinuity waves. These local solutions can 
then be used to replace the Riemann problem as the building block in a 
numerical scheme for a general initial value problem. The advantage of the 
new local solutions is that they permit the use of a much larger step size 
than methods which assume uniform flow in each mesh. 


1.3 Outline of Discussion 

The general initial value problem is treated by decomposing it into 
elementary problems which involve a shock, an expansion fan, and a contact 
line. Chapter 2 is devoted to the problem with one shock wave. A typical 
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case is that in which a piston moves into x > 0 with given initial data for 
x > 0. The piston can later be replaced by a contact line in the initial 
value problem. However, it is found that the system of equations for the 
derivatives obtained by perturbing the shock conditions becomes singular when 
the shock strength approaches zero and, therefore, they are unsuitable for use 
with weak shocks. The system of equations can be regrouped such that a common 
factor which vanishes for zero shock strength is cancelled analytically. For 
shocks of zero strength, the curvature of the shock differs from that of the 
characteristic ahead and behind the shock because of a jump in the curvature 
of the streamline which remains continuous in slope. Solutions are 
constructed by assuming regular power series expansions for regions ahead of 
and behind each discontinuity wave. 

Chapter 3 describes results for problems involving an expansion fan. 

The typical problem is that of a piston at x = 0, t = 0, which moves away 
from the gas lying in x > 0. As before, regular power series are used for 
the regions ahead of and behind the expansion fan, and a special expansion is 
employed for the centered expansion fan region. In this region, the usual 
governing equations for one-dimensional flow are rewritten in Riemann 
invariant form in order to uncouple the variables. 

For an arbitary initial value problem, a shock wave, a centered expansion 
wave, and a contact line will exist simultaneously. The preceding two 
solutions for shock waves and centered expansion waves can be combined by 
regarding the piston path as a contact line. The results of this 
generalization are presented in Chapter 4. 

Elementary solutions are given in Chapter 5 which include flow with 
chemical reactions. The cases of a strong detonation wave, a Chapman-Jouguet 
detonation wave followed by an expansion fan, and a Chapman-Jouguet 
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deflagration with a precompression shock wave are studied. These solutions 
can be combined with the results in Chapters 2 and 3 to describe the complete 
solution for initial value problems including chemical reactions. 

Except for the material in Chapter 5, explicit formulas for the speed and 
acceleration of the discontinuity waves and the spatial and time derivatives 
of the flow properties behind the waves are obtained in each chapter. In the 
cases including chemical reaction, the appropriate systems of equations are 
established, but explicit evaluation of flow properties is omitted for 
brevity. 
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CHAPTER 2 


RELATIONSHIP ACROSS ONE-DIMENSIONAL UNSTEADY SHOCK WAVES 


Figure 1 shows the initial value problem induced by a forward moving 
piston x p * F(t) with F(0) = u p > 0 and F(0) = ap. Initial data which 
have continuous right-sided first derivatives at t = 0 for x > 0 are 
prescribed. The problem here is to determine the flow field for small t 
near x = 0; that is, to find the flow quantities and their first derivatives 
behind the shock, region I, and the shock path. Quantities in region II ahead 
of the shock wave are completely defined by the initial data and can be 
determined independently of the motion of the piston. 



Figure 1.- Shock wave. 
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2.1 Governing Equations 


The differential equations governing one-dimensional, Isentropic flow of 
the medium, except at discontinuities, are as follows: 

p + pu + up =0 conservation of mass (2.1) 

t x x 

p(u,_ + uu„) + p v = 0 conservation of momentum (2.2) 

L X X 

p^ + up x = ~ (p t + up x ) conservation of energy (2.3) 

c 

where p is the density, u is the velocity, p is the pressure, c is the 
speed of sound, and where x, t are independent distance and time variables, 
respectively. Subscripts denote partial derivatives with respect to x 
or t. 


2.2 Shock Conditions 


Across a shock, the jump conditions are given as follows [8]: 


P l v l 


P 2 V 2 


m 


p^ + mvj =* p2 + mV2 

1 2 1 2 

h l + 2 V 1 = h 2 + 2 v 2 


(2.4) 

(2.5) 

( 2 . 6 ) 


where the subscripts 1 or 2 indicate the region behind or ahead of the shock. 
The notation v j , j = 2,1 stands for the velocity of the fluid relative to 
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the shock in front of and behind the shock, respectively; that is, 

V j ” U j ” u s^>* ^ " 1,2. Enthalpy is denoted by h and m is the mass flow 

rate. 

Before perturbing the shock conditions, note that for weak shocks these 
equations become linear equations for u^ - u^ , p^ - p^ , p£ “ Pj. The 
determinant of the system of linear equations vanishes as the shock strength 

goes to zero. Therefore, for weak shocks the solution of these equations will 

be in terms of ratios of small numbers and will Induce inaccuracy in any 
numerical scheme. In order to avoid this difficulty, the basic equations 
should be recombined to cancel out a common factor which vanishes as the shock 
strength goes to zero. With the ideal gas relationship h = yRT/(y ~ 1), 
equation (2.6) becomes 




2yR 
Y - 1 


(Ti 


T 2> 


(2.7) 


Rewriting equation (2.5) using p = pRT and equation (2.4) yields 


v 2 “ V 1 




P2^) “ ( v 2 ” 



where equation (2.7) has been used to substitute for T^. Now the factor 
(V 2 “ Vj) can be cancelled to obtain 


2lfv l v 2 


- (Y - l)(v 0 + 


v l )v 2 


2yRT, 


2c, 
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2 

Finally, the speed of sound relation c 2 = YP 2 /p 2 ~ yRT 2 is u8ec * to arrive at 

(y - l)v 2 2 - (Y + l)v lV2 + 2 c 2 2 = 0 (2.8) 

Equations (2.4), (2.5), and (2.8) are regarded as the basic shock 
conditions. They are, of course, valid for finite shocks, but the determinant 
of the system of linearized equations will now remain finite as the shock 
strength approaches zero. 


2.3 Method of Solution 

The flow quantities for small t behind the shock are related, not only 
to the slope and curvature of the shock and to the initial velocity and 
acceleration of the piston, but also to the flow properties ahead of the 
shock. The following data are provided: 

1. The flow properties for x > 0 at t = 0. Therefore, the first 
partial derivatives with respect to x of the properties are known 
at t = 0. 

2. The velocity u and acceleration a of the piston at t = 0. 

r r 

The piston path x p (t) is then, to second order in t, 

1 9 

x p (t) = u t + j a p t z (2.9) 

In regions I and II, the solutions are regular. Hence, their spatial and 
time derivatives exist as x and t approach 0 + in each region. If Q 
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represents any flow quantity and j - 1,2 represents the regions, then the 
power series expansion about (0,0) of any quantity is 


8 8 

Q (x,t) - Qj (0,0) + x — Qj (0,0) + c Q t — Q j ( 0 , 0 ) 

x c o 

+ higher order terms 


or 


Qj(x,t) = Qj ^ + x Qjx + C 0 fc Q jt + hi 2 her order terms (A) 


, n , 3 t 9 Q 1 

where Qj ■ Qj(0,0), Qj x ” g^Qj(0,0), and Qj fc » — jj*}- at (0,0). The 

speed of sound at (0,0) in region II is taken as the reference speed Cq. Let 

the shock front velocity and acceleration at time t “ 0 be denoted by u 

3 

and a„, respectively. One has 
s 


x g (t) - x g (0)t + “ x 8 (0)t 2 + 0(t 3 ), 


and 


u 8 (t) “ u 8 + a g t + 0(t*) 


(B) 


The procedure used in the following derivation is outlined as follows: 

1. Substitution of the regular power series expansion into the governing 
equations (2.1) to (2.3) and use of the initial data 

u 2^’ p 2^ 0) ’ P 2^’ u 2x* P 2x’ and P 2x’ tD define the flow 
properties in region II, the region ahead of the shock. 

2. Substitution of the expansion series (A) and (B) into the shock 

conditions (2.4), (2.5), and (2.8) and the governing equations (2.1) 
to (2.3), and then comparison of the coefficients of like powers of 
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t to define the velocity and acceleration of the shock and the flow 
properties behind the shock. The boundary condition on the piston 
path which must be satisfied is 

dx 

u i ( v t} = IT' “ u p + V (2 * 10) 

The results obtained for the leading terms should be exactly the same as 
those for a shock separating two uniform regions. The detailed results for 
the acceleration of the shock and the derivatives of the flow properties are 
derived in the next section. 


2.4 Explicit Solutions 

The spatial derivatives of the flow properties for x = 0 + at t = 0, 
that is, and Q2 X t * ie ex P ans i° n series (A), are given. 

Substituting the power series expansion (A) into equations (2.1) to (2.3) 
defines the time derivatives in region II in the form 


_1 r <°> . 

P 2t " ~ [P 2 U 2x + p 2x U 2 ) 
0 


(2.11) 


-1 / p 2x (0) 

U 2 t ~~ r oW^ + “ 2 U2xl 


( 2 . 12 ) 


p* ■ * U 2 (0 , pJ 


(2.13) 


11 



Thus the flow properties and all their x- and t- derivatives in region II 
are defined. 

Now we proceed to the determination of quantities in region I. First, 
the expansion series (A) are substituted into the shock conditions (2.4), 

(2.5) and (2.8), and the results are evaluated on the shock path (B). Then 
the constant terms and the coefficients of t are equated to zero in each 
equation. If it is noted that the relative velocities Vj are 

v 2 ” ( U 2 (0> “ u s) + ( u 2x u s + U 2t C 0 " a s) fc (2 ‘ 14) 

V 1 = ( u l <0) " u s ) + ( u lx u s + u lt c 0 " a s^ (2 * 15) 

then the constant terms yield the standard equations for a straight shock 
front: 


(0) r (0) n (0) r (0) y 

>1 (“l u s) ■ P 2 (“2 -V 


( 2 . 16 ) 


(0) . (0)( (0) w 'i K'JJ \z 

P 1 + p 2 ( U 1 " U s^ U 2 'V = p 2 + P 2 ^ U 2 ~V (2,17) 


( 0 ) 


( 0 ) „ ( 0 ) r „ ( 0 ) 


<Y + iXu / 01 - u s )(u 2 w - u s ) - (y - l)(u 2 w - - 2c 


( 0 ) 


( 0 ) 


( 0 )' 


(2.18) 


The coefficients of t yield three more equations relating the derivatives of 
the quantities ahead and behind the shock front to its initial acceleration: 
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<p lx"s + p lt c O ) ( u l <0> ■ u s) + p l <0)(u lx u s + u lt c O ‘ V 


<0 2x"s + p 2t c 0>( u 2 <0> - + e 2 <0) (u 2x"s + u 2t c 0 ‘ V 


(p lx“s + p lt c O > + p 2 <0) ( u l <0> - u si (u 2x u 8 + u 2t c 0 ‘ a s ) 


+ p 2 (0) (» 2 <0) - u s X“i x “ 8 + u i t c o - V 


+ (p 2x u s + p 2t c 0 > ( u l ( ° > ‘ u a H» 2 <0) - »„) 


< p 2x u s + p 2t c 0 ) 


+ 2p< 0) (u< 0) - u s )( U2x u s + u 2t c Q - a s ) 


+ <p 2x u s + p 2t c 0 ) t u 2 < ° > ‘ 


(Y + l)(u 1 (0) - u s )(u 2 x u s + u 2 t c Q - a s ) 


+ (Y + l)(u 2 <0) - » s )Cu u u s + u 2 t c 0 - » s ) 


- 2(y - l)(u 2 (0> - u s )(u 2x u s + u 2t c 0 - a B ) 


= 2c, 


( 0 )' 


P 2x U s + P 2t C 0 p 2x U $ + p 2t C 0 


(0) 


(0) 


( 2 . 19 ) 


( 2 . 20 ) 


( 2 . 21 ) 


13 



Now the boundary condition, equation (2.10), at t = 0, requires that 


Ui (0) " U r 


( 2 . 22 ) 


Then equation (2.18) defines the velocity of the shock as a function of the 
piston velocity and flow properties at (0,0). This is the well-known result 
for a shock separating two uniform regions. A quadratic equation on u_ is 
obtained in the form 


2u s 2 + [ (Y - 3 )u 2 (°) - (y + 1 ) u p] u s + (Y + l)u p u. 


( 0 ) 


(2.23) 


- ( Y - l)u <°> - 2 c 9 (0)2 = 0 


from which u„ is determined as 
s 


u p ■ JL D u 2 <0> 1 ■‘-P- J ( u 2 <0) - u p, 2 + 


4c 


(0)\2 


+ 1 


(2.24) 


In equation (2.24), the + sign is for a forward facing shock, that is, with 

the gas on the right-hand side, x > 0, and the - sign is for a backward 

facing shock, with the gas on the left-hand side, x < 0. Now using the shock 

2 

speed from equation (2.24) in equations (2.16) and (2.17) with c » Yp/p» the 
leading terms of the flow properties behind the shock can be found. It is 
easily seen that those results are exactly the same as the well-known results 
for a shock separating two uniform regions and are given by: 
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(2.25) 


<0) r (0) x 

... p o f- u 9 “ U J 

(0) _ 2 2 s 


'1 


(u - u ) 
P s 


P,<°) ■ P 2 (0 > + P 2 <°W 0) - u )(u <°> - u) 


(2.26) 


At this stage of the analysis all the leading terms Qj^®) are known. 

To determine the first-order terms Q^ x and Q^ t » the expansion series 
(A) are substituted into equations (2.1) to (2.3). The resulting three 
equations are evaluated in region I using the fact that u^^ = U p 
(eq. 2.22). Similarly, the series (A) is substituted into equation (2.10). 
This yields 


p, u + p, c„ 
lx p It 0 


• -P 


( 0 ) 


u 


lx 


'lx 

u lx u p + u lt c 0 = -Joj 
P 1 

Plx u p + Plt c 0 = -Cl ( ° )2 Pl (0) Uix 


u lx u p + u lt c 0 = a p 


(2.27) 


(2.28) 


(2.29) 


(2.30) 


Equations (2.19) to (2.21) and (2.27) to (2.30) represent a system of seven 
equations for the six derivatives of the flow quantities in region I at the 
point (0,0) and the acceleration of the shock a„. 

o 

To solve these seven equations, it is noted that equations (2.27) to 
(2.30) can be used in equations (2.20) and (2.21) to derive a set of linear 
equations on a g and u^ x as follows: 
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(2.31) 


, <°) ^ u 2 ( 1 ■ u s r (o ) 2 , ,2] 

-fu ' - u ) c, + (u - u ) 

1 2 p J u p - u s L 1 P s' J 


ra n 
s 


2 % (U 2 (0) - U S ) + D 

•“»)*,! 1 l“ z ‘ 01 - «.) (“a 1 ”’ »,) j 


_ U lx^ 


.^ 2 (0) ‘ U s )a p + E . 


The quantities D and E in equation (2.31) are defined as 


0 “ -7olf(P2x u s + p 2t c 0> + P2 (0) < u 2x u s + u 2t c o) [( u 2 (0) " u s ) + ( u 2 (0) " u p )l 
P 2 k 

+ (p 2 x u s + p 2 t c 0 )(u 2 (°) - u s )(u 2 (°) - U p )j 


E ■7Tt( (u 2x u s + u 2t c 0>[ ( l' * 1)(u p - "s> - 2 <1" 1 >t“ 2 (0> - u s )] 


2 c ,<°> 2 


' P 2x U s + P 2t C 0 

p 2x u s + p 2t c 0* 


n (0) 
P 2 

1 

0 

(M 

Q. 

) 


Then the shock curvature a 0 follows as 

8 


a s 3 


(u 2 (0) - u s )[(u p - u s ) 2 - c / 0 ’ 2 ] . (D-E)(u p - « s ) 2 - C, 

«.)« 

■ 


( 0 )‘ 


, (0) . (0) 4 . (0) 

( U 2 - V C 1 ■ 7 TT ( u 2 


/ X 2 ( 0 ) 

(u - u ) + c 

p s 1 


and u^ x is 


Y ♦ 1 P 


u lx 


•„<“» - vh (o1 - " P » - vtt '“i 101 -■.)]* ;,r : ' "- 1 <“ - e ' - r±r »(«. - ■ 

, ( 0 ) , ( 0) 2 4 , ( 0 ) J, .2 ( 0) 2 

( U 2 - U p^ C l ~ T (U 2 -“sTV'V +C 1 


(2.32) 


(2.33) 


(2.34) 


(2.35) 
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The denominator in these expressions is rewritten in the form 


( 0 ) 


I " Y (0) 


( 0 ) 


- < C 1 K*p -V + M (U 2 " U s» + TTT <V ' - u s )( V V > 


Now, Up - u g = and - u g = V£« Therefore, since 

p.v 1 =p 9 v 0 , u -u and u„^^ - u have the same sign, this denominator 
11 4 z p s z s 

cannot vanish. 

Now use of equations (2.27) to (2.30) in equation (2.19) allows p^ to 
be determined as 


(p 2 X U s + P 2 t C 0 ) ( U 2 (0) * U S } + p 2 lU,(u 2x U s + Vo 1 + ^ P / U ' - p 2 W N “ p l' U ' a p + 2p 2 <0) ( u 2 <0> ‘ U S ) U 1 


( 0 ) 


( 0 ) ( 0 ), ( 0 ) 


(u - u )‘ 
P s 


(2.36) 


and the solution for p lx is directly expressed with equations (2.28) and 
(2.30) as 


ei* ■ -vi (0) <2 - 37 > 

Finally, equations (2.27) to (2.29) are used to define the time derivatives of 
the flow properties at (0,0) in the form 


a - u u 
P lx p 


l lt 


It 


- (p l (0)u lx + p lxV 
c 0 


(2.38) 


(2.39) 
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Pit = 


(2.40) 



The shock problem with a piston moving into the flow has now been solved 
completely for small t, near x » 0. First, equations (2.11) to (2.13) give 
the time derivatives for the region ahead of the shocks, which are 
completely defined by the initial data, and <?2x> independent of the 

piston motion. Next, equation (2.24) gives the velocity of the shock, and 
then equations (2.22), (2.25), and (2.26) define the leading terms of the flow 
properties behind the shock; they are exactly the same as the well-known 
results for a shock separating two uniform regions. Equation (2.34) yields 
the acceleration of the shock. The derivatives of the flow properties behind 
the shock at the point (0,0), that is, u lx , p lx , p lx , u lt , p lt , p lfc , are 
determined in equations (2.35) to (2.40), respectively. Thus, the velocity 
and acceleration of the shock as well as the first three terms in the power 
series expansion of all the flow quantities behind the shock have been defined 
in this chapter. 
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CHAPTER 3 


CENTERED EXPANSION WAVES 


The typical problem involving a centered expansion wave is that in which 
a piston starts to move at t = 0 away from a gas which lies in x > 0. The 
problem in this investigation is to determine the relation of the initial flow 
properties in x > 0 and the velocity and acceleration of the piston to the 
flow properties in the expansion fan, the boundary of the fan, and the flow 
properties behind the fan. 


3.1 Governing Equations 

The governing equations for one-dimensional flow are equations (2.1) to 

(2.3). These equations can be rewritten in order to uncouple the variables. 

By using Riemann's invariant approach, the characteristics of the differential 

equations of one-dimensional isentropic flow can be found. Combining equations 

2 

(2.1), (2.3) and c =■ yp/p, yields 

Y pu x + p t + up x " 0 
Using equation (2.2) one has 

c ( puu + pu + p ) =* 0 

X LX 
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Adding and subtracting these two equations yields 


p t + (u + c)p x + pc[u fc + (u + c)u x J = 0 (3.1) 

p t + (u - c)p x - pc[u t + (u - c)u x ] - 0 (3.2) 


2 

Total differentiation of both sides of c = Yp/p with respect to t and use 
of equation (2.3) gives 


dc _ Y “ 1 dp 

dt p dt 


(3.3) 


This relationship can be applied in equations (3.1) and (3.2), to rewrite the 
system as 

(ft + (u + c) ir)( u + 7nr c ) 

(|_ + (u . c) |_)( u . _i_ c ) 

Equation (3.3), with the relationship c i = yp/p> can be written as 



Now equations (3.4) and (3.5) express the relationship between u and c, and 
equation (3.6) can be used to define the pressure p. These three equations 
can then be used to replace the original system of equations (2.1) to (2.3). 


= 0 (3.4) 

= 0 (3.5) 
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3.2 Method of Solution 

Define region I as the region behind the expansion fan. Region II is the 
region ahead of the expansion fan where the flow properties are not affected 
by the motion of the piston. The expansion fan itself is called region III 
(Fig. 2). 



Figure 2.- Expansion fan. 

OA is the characteristic line that separates regions III and II, and OB is the 
characteristic line that separates regions III and I. For a forward facing 
expansion wave, OA and OB are forward characteristic lines C + , on which dx/dt 
= u + c. For a backward-facing expansion wave, they are C - characteristic 
lines with dx/dt = u - c. 

Near the center of the expansion fan, the flow quantities must jump from 
their values at OA to those at OB. Locally, they behave as a simple expansion 
wave and are functions of the slope of the characteristic lines at the center 
(in the x, cgt plane). It is convenient to introduce the angular coordinate 
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n ” x / c q t as an independent variable to replace x. With n and t as 
independent variables in region III, any flow quantity, Qg(x,t), can be 
expressed as a power series in t with coefficients which are functions of 
ns 


Q 3 (x,t) = Q 3 (0) (-^ r ) + + hl g her order terms 

«* Q 3 ^ (n) + Q 3 (1) (n) c^t + higher order terms (C) 

A characteristic curve will be specified by a parameter £, where 5 is the 
slope of the curve at the origin (see fig. 3). 



Figure 3.- Geometry of the expansion. 


Thus, the family of characteristics is given by x * x c (£,t) where 

x (5,t) = £c_t + k(5)(c.t) 2 + 0(t 3 ) in which k(£) => d 2 x /2d(c n t) 2 
c u u c u 


t=0 


is 


to be determined in the course of the analysis. If ^ anc * ^ are 

slopes of the first (OA) and last (OB) characteristics at the origin, then the 

family of characteristics is given by < K < Ky* Note that, by definition. 
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the angular coordinate is related to the parameter 5 on a characteristic 
through 


x c (S,t) 

n = — = 5 + k(£)c n t + higher order terms 

c Q t U 


(3.7) 


3.3 Outline of the Steps 

The mathematical procedure is to substitute expansions (A), (B) and (C) 
into the governing equations and then to compare the coefficients of like 
powers of t. The outline of the steps follows: 

1. Determination the flow properties in region II by substituting series 
(A) into equations (2.1) to (2.3). This will then determine the 
first characteristic line in region III, OA, and its slope and 
curvature, by use of the boundary condition which requires that the 
flow properties are continuous across OA. 

2. Determination of the flow properties in the expansion fan, 
region III, by substitution of series (C) into the governing 
equations (3.4) to (3.6). 

3. Determination of the last characteristic line in III, OB, and its 
slope and curvature. 

4. Application of the continuity condition along OB and determination of 
the first derivatives of the flow properties in region I, which must 
be compatible with the velocity and acceleration of the piston. 

The results for centered expansion waves, the velocity and acceleration 
of the waves, and the gradients of the flow properties are derived in next 
section. 
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3.4 Explicit Solutions 


The time derivatives for flow properties in region II are obtained the 
same way as for the shock problem in Chapter II. The equation of the 
characteristic OA is 


x 


= x c (C 2 ,t) - 5 2 c Q t + k(? 2 )(c Q t) 2 + ... 


(3.7a) 


The flow properties are required to be continuous across OA so that 


Q 


2 


+ Q 2x X c (5 2’ t) + Q 2tV = Q 3 ( x c ^ 2 ,t),t) 
= Q 3 (0) (n 2 ) + q 3 (1) (n 2 )c Q t 


2 

where n 2 = x c ^2 ,t ^' C 0 t = ? 2 + k(5 2 >c Q t + 0(t ) from equation (3.7). Thus, 
using Taylor series expansion, this equation becomes (up through linear terms 
in t) , 


( 0 ) 


+ QfcihV + Q 2t c 0 t 


q 3 (0) (5 2 ) + k( 5 2 )c 0 t q 3 (0, '<e 2 ) + Q 3 (1) <? 2 >V 


(3.8) 


Now, comparing the coefficients of like powers of t in equation (3.8), it is 
found that 




(3.9) 
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and 


, (1) <5 2 > -<W ; 2 + «2t -k(5 2 )Q 3 (0) '<5 2 ) 


(3.10) 


Equations (3.9) and (3.10) represent boundary conditions on the functions 
Q < 0 >(„) and (h) at the point n = £ 2 and will be used in the 

following steps. 

By the definition of the characteristic lines, dx/dt = u + c for 
forward-facing expansion waves, and dx/dt = u - c for backward-facing 
expansion waves. Thus 


^2 C 0 + 2k ^2 )c 0 t = u 2 + u 2x ^2 c 0 t + u 2t C 0 t 


* (C 2 (0) + c 2x 5 2 c 0 t + c 2t c 0 t) 


(3.11) 


Comparing coefficients of the same power of t, one obtains 


u (0) ± c (0) 
2 1 2 


'0 


and 


(3.12) 


, rr U 2x * C 2x r U 2t ± C 2t 
k(S 2 ) = 5 2 + 


'0 


'0 


(3.13) 


Hence, the slope and curvature of the head characteristic OA are defined. 

To proceed with the second step in the analysis, the special expansion 
(C) is used in region III. The spatial and time derivatives are as follows: 

Q 3t (n,t) " -Q 3 (0r < n > t ~ Q 3 (1>,(n)nc 0 + Q 3 (1>(n)c 0 (3 * U) 

and 
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Q 3x (n,t) “ Q 3 ( ° ),(n) TT + Q 3 (1),(n) 0.15) 

Substituting these equations into the governing equations (3.4) and (3.5) and 
then comparing the coefficients of the same power of t, one gets 

j^ _c 0 T ' + u 3^^( T l) + c 3 ^(h)J£u 3 ^) (r)) + - j ^ (q)J = 0 (3.16) 

|^ _c 0'n + u 3 ^°^ ( n) - c 3^(ti) u 3 (°) (n) - ~~T = 0 

^- c 0 t) + u 3 (°)(n) - c 3 (°)(n)J £u 3 ( l ) (n) - — ^ c 3 ^ 1 ^ ' ( n )J + u 3 ( 1 ) ( n) Jjc 0 + u^^'U) 

‘ 7 TT c 3 ( ° ) '<’i)j ’ C 3 ( 1 ) <n)[^ r |- T c 0 + u 3 <°>'(n) - — - c 3 (°)'(„)J - o 

1 . 

[■■ C 01 + u 3 (°)(ti) + c 3 (°)(T,,][u 3 (l)‘( n ) + — c 3 < 1 > * ( n )J + u 3 < l )(i,)[c 0 + u 3 < 0 >* (ti) 

+ ~T c 3 (0) '<i)] + C 3 (1) <n)|~r-T c o + u 3 (0) 'U) + ~ T c 3 <0) * (n)] * o 

From equations (3.16) and (3.17), two sets of results are obtained: 

(3.20) 

(3.21) 



(3.17) 


(3.18) 


(3.19) 
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in which the upper signs correspond to forward-facing expansion waves and the 
lower signs correspond to backward-facing expansion waves. 

Integration of equation (3.21) and application of boundary condition 
(3.9) at n =* (>2 yields 


u.^^(n) + — — r Co(°)(n) = const 
j Y - 1 J 




(3.22) 


If a constant B is defined by 

(0) (0) 

u c 

2 2 2 

B = + + 

c Y - 1 c 

0 0 

then the terms c^^(n) and u^^(n) can be written, using 
equations (3.20) and (3.22) as 


(3.23) 


co^(n) = ± (n ± B)cq (3.24) 

J Y + i u 

u 3 (0) (^) = ttt( T1 b ) c o (3,25) 

Equations (3.20) and (3.22) are now substituted into equations (3.18) and 
(3.19), which yields 

u 3 (1) (n)[co ± y-fr c 3 (0),(t i ) ] * ~TJ C 3 (1 ^( T ' ) [ C 0 1 2c 3^°^ ' ( n)J =0 < 3 * 26 > 

+ 2 c 3 (°)(n) Ju 3 ^^ (ri) + i c 3^ (h)J + J^ u 3^ ^ ( h) + ^ j f c 3^ (h)j c o = 0 

(3.27) 


Substitution of equations (3.24) and (3.25) into equation (3.26) results in 
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u-j^Cn) = ± Acg^^Cn) 


(3.28) 


where 


A = ~2(3y ~ 1) 

(Y “ D(Y + 5) 

It follows that 

u 3 ( 1 ), (n) = ± Ac 3 ( 1 ), (n) 

which can then be substituted into equation (3.27) to obtain 
2c 3 ^ 0 ^(ri)c3^ 1 ^ , (n) + c 3 ^ 1 ^(n)c 0 - 0 


(3.29) 


or 


c (D’ (n) (D (n) . o 

c 3 2 (y - 1) n ± B 3 ^ n; 

The general solution for c 3 ^(n) is given by 

Y fl 


c 3 ^^(ti) = c- 




2(y-D 


(3.30) 


where, using the boundary condition (3.10), 

c 3^ ^^2^ = ^2 c 2x + c 2t + y + 1 k ^2^ c 0 (3.31) 

and then u 3 ^^(n) is known from equation (3.28). Of course, from equation 
(3.28) it follows that u^^^E^) “ ± Ac^^CE^), an ^ this must be consistent 
with the boundary condition (3.10) on u^^^E^) an< * c 3^^^2^* This has 
been proved by the author by applying the characteristic equations (3.4) and 
(3.5) in region II and using equation (3.12) for This yields a 
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relationship between ^2 n 2x + U 2t anc * ^2 C 2x + C 2t* ^ rom it follows 

that * ± Ac 3^^ (^ 2 ) • 

Now equation (3.6) represents the relationship between u, c, and p. 
Since u and c are already known, this equation can be used to calculate 
pressure. Substituting the expansion series (C) into equation (3.6) and 
comparing the coefficients of like powers of t yields two equations. The 
first one is 

(3.32) 


- 1 


c 3^(t))P 3^ (b) 




= 2p 3 ^ ( ti)c 3 ^ 0 ) (q) 


U3 (0| U> 


and the second one is 

(n)|jy c 3 (°) ( Ti)J + P 3 ^ 1 ^ ( rj) j~±2c3^ 0 ^ ( 11)03^°^ ( ri) + - C3^^(q)c Q j 

+ P 3 ^(n)j* ^ ~ - c 3 ^ 1 ^ ( n)c 3 t°^ ( 13) + U3^)(r))c3(°)(q)j 

- 2p 3 <°){t 1 )[tc 3< 1 ) ‘(T,)c 3 ( 0 )( n ) + c 3 U)(q)c 0 + u 3 ( 1 )(T 1 )c 3 (0) , ( n )] - 0 (3.33) 


Since the factor -ri + u 


( 0 ) 


(n)/c, 


+ c 3 < 0) /c 0 is not zero, it can be 


\ 1 I / / V-Q 

cancelled in equation (3.32) resulting in a first-order ordinary differential 

( 0 ), 


equation for p^ (n): 


P 3 ( 0 ) '(n) - 2Y 


Y - 1 n ± b 


p 3 ^°^(n) =■ 0 


Integrating this differential equation and using the boundary condition (3.9), 
the solution for P 3 ^^(n) can be obtained as 
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P3^U) * P2^| 


A differential equation for ^ (n) can also be obtained by rewriting 

equation (3.33) with the aid of equations (3.24), (3.28), (3.30) and (3.34) 
in the form: 




The general solution of this inhomogeneous equation is 




3( y+1 ) 

2 ( y- 1 ) 


P3 », (n) . ^ ^ 0) “iMI * «- ‘ '» rl 


U 2 t B) 


(5y+l) 

2(y-l) 


The integration constant, G, can be written as 

u ^(f ) 

n (1)/- \ 1 . y(y + l)(y + 5) 3 (o) 1 

G ■ Pa <e 2 ) 3^i (t . i )(3y . i) c Pa TT 

(«, t „r> <«* * 


so that solution p^ '(n) is represented by 


{ u (1) fr \ D {a) -3y*i 

t(t ♦ l)(r ♦ 5) u 3 P 2 /„t B\2 ( t -1) 

(V- 1)(3t - lj c Q ^ * B t bJ 


+ 1 + p 3 (1)(5 2 ) 


I'e * B J 


and p^ '(5 ^ ) is defined by the boundary condition (3.10) in terms of known 


quantities as 
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( 1 ) 


(c 2 ) 


P 2x 5 2 + P 2t " k( ^2 ) p 3 


(0) 


(? 2 ) 


Proceeding further, the continuity equation (2.1) can be used to define 
the density. The procedure is the same as before: substituting expansion 

series (G) and its derivatives with respect to x and t into equation (2.1) 
and then comparing coefficients of like powers of t. Two equations are 
obtained: 


P3^ (t)) | 


“ 3 (0) ^ 

c o J 


. u ( 0 ) '(t,) 

+ p 3 0 (t)) 

0 


= 0 


(3.36) 


and 


P3^ (t))[-t)Co + U3^U)] + P3^^(n)[ c o + u 3^ U)] + 03 ^) (r))p3^ 0 )(n) 


+ U3^)(q)p3^^ (t)) * 0 


(3.37) 


Now, solving for p^^(n) from equation (3.36), with the help of equation 
(3.25) and the boundary condition (3.9), one obtains 

2 

P3 (0) U) ■ P2 (0) (~ J^ 1 " 1 <3.38> 

Similarily, solving equation (3.37) using equations (3.25), (3.28), (3.30), 
and (3.38), the solution of p^^(n) takes the form: 
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P3^(*l) 


(y + 1 )(y + s ) 


o 2 <°> 


(y - l)(-3y + 1) 


l z t B 




-3y+l 

2(y-l) 



(3.39) 


where known from the boundary condition (3.10). 

At this stage, all the flow properties in the fan region III are defined. 
The next step is to define the tail of the fan OB; that is, its slope and 
curvature at t m 0. Since OB is also a characteristic line, with dx/dt 
* u ± c on OB, the flow properties must be continuous across it. The 
boundary conditions on OB are similar to those on OA. First denote the slope 
of the OB at origin by 5^. The characteristic line OB can be written as 
follows: 

x - ^(.(Sj.t) - CjCQt + k(5 1 )(c Q t) 2 (3.40) 

Now, the boundary condition along OB is expressed as before as 

Qi (0) + Qix^ot + Qit^ 

- Q 3 (0) (5 1 ) + k(5 1 )c 0 tQ 3 (°>'(5 1 ) + Q 3 (1) (5 1 )c 0 t 

By comparing the coefficients of the same power of t, two equations are 
derived: 




(3.41) 


*lx 5 l + Q lt 


Q 3 (1) (5 1 ) + k(5 1 )Q 3 ( ° r (5 1 ) 


(3.42) 
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The first of these can be used to define the slope of OB; that is, if Q is 
identified as u, than equations (3.41) and (3.25) define using the fact 

that ^ (see eq. 2.22): 




, u 

Y + 1 p Y “ 1 

2 c o 2 


B 


(3.43) 


The constant B is defined in equation (3.23). Substituting into equation 

(3.20) and appling equation (3.41) yields 


(0), p n x <°)/ r \ (0) . (0) 

U 3 (5 1> * c 3 'V "l * °1 


0 


'0 


(3.43a) 


from which 


Cl ( °) “ * -- -- -- (u p ± Bc 0 ) 


(3.44) 


By definition, on the characteristic line dx/dt - u db c. This implies that 


^Cq + 2k(5 1 )CQ 2 t - U 3 (°)(Sl) + k(5 1 )c Q t U 3 C0) , (5l ) + ^^(Sj.JCgt 
± (c 3 (°>( ?1 ) + k(? 1 )c 0 t + c 3 ^ 1 ^(? 1 )c 0 t) 


With equations (3.24), (3.25), (3.28), (3.30) and (3.43), every term in this 
equation is known except k(£^). Therefore 


k(^) 


U) f , ) 

(Y + 1 ) ( Y - 3) U 3 


2(3y - 1) 


(3.45) 


'0 


where u^^(S^) is defined in equation (3.28). This completes the third 
step in the analysis as outlined previously. 
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The final step consists in using the slope and curvature of OB and the 
boundary conditions on the piston path to define the flow properties in 
region I. First, substituting the expansion series (A) into the governing 
equations (3.4) and (3.5) yields 


u it c o + ( u i (0) + c i (0) ) u ix + c it c o + K (0) + c i (0) ) c ix " 0 < 3 - 46 > 


and 


u lt c 0 + ( U 1 


( 0 ) _ 



c lt c 0 




0 (3.47) 


Adding and subtracting equation (3.47) from (3.46), this system of equations 
becomes 


u lt c 0 + u i^ u i~ + r c/°)ci„ » 0 


l lx 


Y - 1 


lx 


u lx c l (0) + 7TT c lt c 0 + TTT u l (0)c lx “ 0 


(3.48) 

(3.49) 


along with these, the boundary conditions on OB , equation (3.42), imply that 
u lt + u lx C l “ k <5 1 )u 3 (0), (? 1 ) + u 3 (1) (5 1 ) (3.50) 

C lt + C lx 5 l “ k(5 l )c 3 <0) ' (5 l ) + c 3 (1)(? l ) (3,51) 


However, it is easily seen using equation (3.43a) that the determinant of the 
system (3.48) to (3.51) vanishes. A fourth independent relation can be found 
in equation (2.30): u i x u p + u it c 0 “ a p* Using this yields 
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(3.52) 


c lx 


+ 


a 

P 


u ± Be 
P o 


“lx - |l 


4(Y+1>C ° uj<l> (ei )- 2 


( 3 y - 1 )(y - 1) 


- 1 P 


1 


I u ± Be 

p 0 


and 


c lt = + 


|( r + l)(r. - l) u (D( e ) . jh 

1 «®r- 3 1 VVl 


J lt 


4 ( y + l)u p 
I ( 3 Y - i)(y - 1) 


u 3 (1) <5!> -ttt Vijrr 


Be. 


Substituting expansion series (A) into equation (3.6) for region I, 
yields 

(0) 

2y P 1 


Plt c 0 + Plx u p “ T^T “(O)" (c lt c 0 + c lx u p> 


From the boundary condition on OB 


s. ± B\y-1 


pi 


(°) = p 3 (0)( 5l ) = p 2 (0)^_I__j 


and 


(3.53) 


(3.54) 


(3.55) 


(3.56) 


(3.57) 
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"it + p lx 5 l 


(3.58) 


- p 3 (0) '(5 1 )k(5 1 ) + P 3 (1) (C 1 ) 


Therefore, equations (3.56) and (3.58) with equations (3.43) and (3.49) yield 
for p^ x the expression 


YP 


( 0 ) 


( 0 )' 


( 0 ), 


u lx + P 3 V “' (5 l )k(5 l )c O + p 3 U 1 )C 0 


lx 


Y - 1 


(Up ± Bc 0 ) 


(3.59) 


Substitution of equation (3.59) back into equation (3.56) then defines P]^* 
The procedure that defines the density p in this region is the same as 
that which defines the pressure. The equations used are the continuity 
equation and the boundary condition, equations (2.1) and (3.42), with Q 
identified as p. They are as follows: 


Pi _c n + p. u 
It 0 lx p 


„ ( 0 ) 

P 1 u lx 


(3.60) 


p lx 5 l +p lt- p 3 < ° ,,( 5l )k(5 l> + ‘ > 3 <1) < 5 l ) 


(3.61) 


Then the spatial derivative of p can be determined as 


( 0 ) 


U, + p,"" (ej^S^CQ + p, (1) (5,)c, 


lx 


( 0 )' 

'3 U 1 


’ 1^0 


lx 


L i A ( % ± Bc 0> 


(3.62) 


Substitution of equation (3.62) back into equation (3.60) then defines P^ t * 
The problem involving a centered expansion wave, that is, in which a 
piston moves at t « 0 away from the gas, has now been completely solved. 
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The governing equations inside the fan are expressed in Rlemann invariant 
form, equations (3.4) and (3.5), in order to uncouple the variables. This 
form relates u and c, and then the original differential equations are used 
to define p and p. The regular power series (A and B) are used for the 
regions ahead of and behind the expansion fan, and a special expansion (C) is 
employed for the centered expansion fan region. 

The flow properties in the region ahead of the fan depend only upon the 
initial data of the flow and can be expressed in the same manner as for the 
region ahead of a shock. The head characteristic, OA, its slope and 
curvature, are defined by using the boundary conditions, equation (3.8), that 
require the flow properties to be continuous across OA; they are given in 
equations (3.12) and (3.13). The solutions for the flow properties inside the 
fan are in two parts. The leading terms, c^^(n), u^^(n), p^^^n) 
and p^^(n), are defined in equations (3.24), (3.25), (3.34), and (3.38), 
respectively. These are the same as the solutions for a centered expansion 
wave with uniform initial data, and they are functions only of the flow 
properties ahead of the fan. The first order terms c^ * (n)» u^ .^'(n), 
P 3 ^(n), an< * P-j vh), are defined in equations (3.30), (3.28), (3.35) 
and (3.39), respectively. They are also related only to the flow properties 
ahead of the fan. The piston velocity at t ■ 0 first comes into the picture 
when defining the slope of the tail characteristic in equation (3.43). The 
curvature of the tail characteristic, equation (3.45), can be expressed after 
its slope is defined. It is interesting to note that the acceleration of the 
piston does not affect the solution anywhere ahead of and within the expansion 
fan. The leading terms in the region behind the fan are defined by equation 
(3.41), with Q identified as u, c, p, or p. Initial spatial and time 
derivatives of c and u in region I are defined by equations (3.52) to 
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(3.55), and the initial spatial and time derivatives of p and p are 
defined by equations (3.59), (3.62), (3.56) and (3.60). Thus, the various 
series expansions of the flow properties in the fan and behind the fan have 
all been determined. 
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CHAPTER 4 


INITIAL VALUE PROBLEM 


The initial value problem with discontinuous initial data joining two 
uniform regions was studied by Riemann [9] and was used as the canonical 
solution, or the building block, by Chorin [3] for the numerical solution of a 
general initial value problem. In this section, an improved canonical 
solution for the propagation of discontinuous initial data joining two 
nonuniform regions is constructed by including the first order effects of the 
nonuniformities . 

For such an initial value problem, a shock wave, a centered expansion 
wave, and a contact line all exist simultaneously. An expansion wave 
propagates into the higher pressure side, a shock wave propagates into the 
lower pressure side. 

A contact line is a particle path separating the region behind the shock 
from the region behind the expansion wave. Across it the temperatures and the 
densities may be different, but it is necessary that the pressure and fluid 
velocity be the same. The problem is to relate the velocity and acceleration 
of each discontinuity wave (that is, the shock wave, the expansion fan, and 
the contact line), to the initial flow properties and their gradients 
(see fig. 4). 
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c O t Particle 



Figure 4.- Shock wave, expansion fan, and contact line. 


The (~) indicates the flow properties on the higher pressure side, that is, 
the expansion wave side. 


4.1 Method of Solution 

The preceding two solutions for shock waves and for centered expansion 
waves are to be combined by considering the piston path as a contact line, so 
that the fluid velocity in the region behind the shock and in the region 
behind the expansion fan are the same along this line. The condition for 
continuity of the pressure across the contact line gives two equations, which 
define Up and ap, the speed and acceleration of the contact line. 
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4.2 Explicit Formulas 


The condition for continuity of the pressure across the contact line 

gives 


Pl (0) + ( Plx u p + Plt c O> t = Pl (0) + ( Plx u p + Plt c 0> t C 4 * 1 ) 

The constant term and the coefficient of t yield 

Pl (0) . p x (0) (4.2) 

and 

Plx u p + Plt c 0 “ Plx u p + Plt c 0 < 4 ‘ 3 > 


By substituting equation (2.24), the solution for u g , into equation (2.26), 
pj^®) is related to u p from the shock solutions 


Pi 


(0) . p 2 (0) t 


+ l 
4 


P2 


( 0 ) 


[(u 2 (0) - Up) * J 


( u 2 (0) - u p) 2 + (7TT C 2 (0) ) 2 (U 2 (0) ' U P } < 4>4 > 


where the upper sign applies for a forward-facing shock and the lower sign for 
backward-facing shock. 

Similarly, by substituting equations (3.12) and (3.43), the slopes of the 
head and tail characteristics of the expansion fan, into the definition of 
Pj^^ given by equation (3.34), and by the use of (3.23), the definition 
of B, results in a relation between p^^ and u p in the form 
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(4.5) 


Pi 


(o) = ? 2 (o) 


Y - 1 
2 


u + Be 
P 0 

c- 2 (0) 


£r_ 

Y-l 


where the upper sign applies to a backward-facing expansion fan, and the lower 
sign to a forward-facing expansion fan. Equating the expression in equations 
(4.4) and (4.5) , results in a relation which determines UpS 


P2 


(0) 



u p) * J( u 2 




,( 0 ) . 


u d) 




u d) 



. u ~ < 0 >\n 


( 0 ) 


(4.6) 


Equation (4.6) must be solved iteratively for Up. Note that the upper sign 
holds for the case of a right shock/left expansion fan; that is, with the 
pressure in x < 0 higher than that in x > 0 at t = 0. The lower sign is 
for the case of a left shock/right expansion fan; that is, for the higher 
pressure in x > 0. The (~) indicates flow properties corresponding to the 
higher pressure side. After Up has been calculated, all the leading terms 
in the expansions of the flow properties on either side of the contact line 
can be determined by using the equations in Chapters 2 and 3. 

The curvature of the contact line can be found using equation (4.3). 
First, repeating equation (2.29), in the region behind the shock 
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where can be determined by equaitons (3.25), (3.26), and with the 

relation - Yp/p» and u lx are defined by equations (2.25) and 

(2.35). Similarly in the region behind the expansion fan, 


Plx u p + Plt c 0 = " 


(4.8) 


In consequence of equation (4.3), the right hand sides of equations (4.7) and 
(4.8) must be equal. Substitution of the expressions (2.35) for u^ x and 
(3.53) for u^ x then gives 
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This completes the required analysis of this chapter. The determination 
of the contact line has been accomplished by joining the solutions for the 
shock problem and the expansion wave problem together using the conditions of 
continuity of pressure and velocity across the contact line. The initial 
velocity and acceleration of the contact line are given by equations (4.6) and 
(4.9). Thus, the solution of the initial value problem with discontinuous 
initial data joining two uniform regions has been extended to the nonuniform 
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case. The canonical solution for discontinuous initial data joining two 
uniform regions has been used as the basic building block for the random 
choice method and for front tracking. With the results of the present work, 
it is now possible to explicitly include information about nonuniform initial 
data in such a building block. This leads to a numerical procedure which can 
utilize a much larger step size in the vicinity of discontinuities in the 
initial data. The actual implementation of the current results into an 
improved numerical scheme for general initial data is beyond the scope of the 
present paper. 
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CHAPTER 5 


UNSTEADY FLOW PROBLEMS WITH CHEMICAL REACTION 


Flow problems with chemical reactions are discussed in this chapter. The 
problem is to relate the initial flow properties and their gradients to the 
speed and acceleration of the reaction front and other discontinuities. The 
piston moving into the combustible gas for x > 0 at t = 0 is treated 
first. These solutions can then be combined with the shock or expansion fan 
solutions in the preceding chapters using exactly the same procedures as In 
chapter 4 for the initial value problem without chemical reaction. The 
algebra involved in this combination is extensive for the initial value 
problem with chemical reaction; thus, it will not be carried out in the 
present work, and only the separate elementary problems will be treated. In 
addition, the elementary problems themselves will not be solved explicitly. 

For brevity, the solutions will be considered complete where an appropriate 
set of algebraic equations governing the unknown quantities is established. 

The shock conditions across a reaction front are [10]. 

P 1 V 1 “ P 2 V 2 " m (5.1) 

P 1 + mv l = p 2 + mv 2 (5.2) 

h l + 2 V l 2 + q l = h 2 + 2 V 2 2 + q 2 (5.3) 
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where p is the density, p is the pressure of the gas, v is the relative 
velocity of the flow with respect to the shock, h is the enthalpy, 

(h a - ■ - T -— £■ is assumed) and q is the energy formation which can be 
released through chemical reaction. The subscript 1 refers to burned gas and 
the subscript 2 refers to unburned gas (i.e., gas which has not yet undergone 
the chemical reaction). In the present section it is assumed that part of q 
is released instantaneously in an infinitely thin reaction zone and the 
unburned gas is on the right. For the sake of simplicity, it is assumed that 
Yj 53 Y 2 = (The case y^ j 4 Y 2 I s more difficult only because of 
additional algebra). When y^ = Y 2 = Y the reaction can be exothermic (i.e., 
can release energy) only if q 2 > qp that is, for the same pressure and 
density, the total energy and enthalpy of the unburned gas is always greater 
than that of the burned gas. In this model, viscous effects, heat conduction, 
and radiative heat transfer are neglected. 

5.1 Chapman- Jouguet Waves 

For completeness, a brief summary is given in this section of the 
elementary theory of one-dimensional detonation and deflagration waves, and a 
derivation of relations between the variables on the two sides of such waves 
is carried out for later use. The material presented here is essentially as 
found in references 9 to 11. Chapman-Jouguet waves have a special 
significance in many systems. It is therefore of interest to investigate the 
properties of these waves first. With the aid of equations (5.1) and (5.2), 
equation (5.3) can be written as 
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a Hugoniot curve (Fig. 5) 




Figure 5.- Hugoniot curve for detonations and deflagrations. 

The intersection between the Hugoniot curve and the straight line through the 
point (1,1) establishes the final state of the system. Each Hugoniot curve is 
therefore divided into two disconnected branches, an upper branch called the 
detonation branch, and a lower branch called the deflagration branch. 
Combustion waves are termed detonation waves or deflagration waves according 
to the branch of the Hugoniot curve upon which the final condition falls. In 
passing through a detonation wave, the gas is slowed down and its pressure and 
density increase; conversely, in going through a deflagration the gas speeds 
up and expands, and its pressure decreases. 

There are two points on the Hugoniot curve at which the straight line 
from (1,1) is tangent to the curve (points D and C on Fig. 5). These are 
referred to as the Chapman -Jouguet points. The slope of the Hugoniot curves 
is given by 
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(5.6) 
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while a straight line through the point (1,1) with this slope has the equation 
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or 
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(5.8) 


Simultaneous solution of equations (5.5) and (5.8) yields the ordinates of the 
Chapman-Jouguet points as 


it 

± 


1 + o(y 



(5.9) 


and the abscissas of the points as 



(5.10) 


in which the upper signs correspond to the detonation branch and lower signs 
correspond to the deflagration branch. A Chapman-Jouguet wave is defined as 
one for which the pressure and density ratio ahead (region II) and behind 
(region I) are given by equations (5.9) and (5.10). 
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Using the continuity and momentum equations and the relation 

2 

c « YP/p* it can be shown [10] that the Mach number of the flow ahead of the 
Chapman-Jouguet front (region II) relative to the front is given by 





(5.11) 


According to equation (5.11), the initial Mach number always exceeds unity for 
Chapman-Jouguet detonations and lies between zero and unity for Chapman- 
Jouguet deflagrations. In addition [9], it is known that Chapman-Jouguet 
waves correspond to the minimum possible propagation speed for detonations and 
the maximum possible propagation speed for deflagrations; therefore, all 
detonations propagate at supersonic velocities, and all deflagrations 
propagate at subsonic velocities. 


5.2 Explicit Formulas for Chapman-Jouguet Waves 

If the flow fields separated by the Chapman-Jouguet wave are nonuniform, 
the wave will not propagate with a constant speed. To relate the Initial 
acceleration of the wave to the gradients of the initial data, we introduce 
again the power series expansion (A) for flow properties Q near the wave front 
which is taken at x ■ 0 for t * 0; 

Qj(x,t) « Qj^ + xQj x + c o t Qjt + higher order terms (A) 
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Let and a^ be the reaction front speed and acceleration at the origin, 

respectively. The reaction front path, x^, is then defined as 


x d (t) 


u d l + 2 a d fc2 + ° (t3) 


(5.12) 


Substituting (A) into equations (5.9) to (5.11), evaluating along the reaction 
front, and then comparing the coefficients of like powers of t, the leading 
terms of the flow properties behind the reaction front are all defined 
explicitly: 




+ 


qy(y - 1) 



± 



(5.13) 
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(5.14) 




(5.15) 


and the speed of the reaction front is found as 
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In equation (5.16) and elsewhere, the notation (±) is used to indicate the 
direction of propagation of the reaction front. The upper sign corresponds to 
a right-facing front, that is, with the unburned gas on the right, and the 
lower sign corresponds to a left-facing front. The upper and lower signs 
which appear without enclosing parentheses refer to the Chapman- Jouguet points 
for detonation and deflagration, respectively. 

The coefficients of the linear terms in t from the foregoing 
substitution lead to a set of linear equations which couple the derivatives of 
the flow quantities; they are 

1 ± 
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and 


> d = u 2 (°)u 2 (l) ( ^ |c 2 (0) 



(p 2 (D - p 2 (') ) (5.19) 


! ♦ Hill : i> 


In equations (5.17) to (5.19), the following definitions have been introduced: 
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Now it is known that a Chapman- Jouguet reaction front, when observed from the 
burned gas behind it, moves with speed of sound [9,11]. This is the famous 
conclusion made by Jouguet in 1905. This relation gives 


(0) (0) 
u i ■ U d <+) V 


(5.20) 
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and 


U l^ ~ ^7oy( a <l(*) \ Cl {0) Cpi (1) “ (5.21) 

Since and c^^ are already defined, the leading term of the flow 
velocity behind the Chapman- Jouguet reaction front uj^) is determined. Note 
that the region behind the reaction front (Fig. 6) is a regular region similar 
to the region behind the shock treated in Chapter 2. 



Figure 6.- Wave pattern of Chapman-Jouguet detonation front. 

Suppose the Chapman-Jouguet wave is initiated by a piston motion, with piston 

•ff ^ 

initial velocity and acceleration u and a , respectively. The piston 

r r 

path is then x p = u p t + 1/2 a p t . The boundary condition on the piston 
path, equation (2.10), requires that 
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(5.22) 


* * 

a = u, u + u-.c A 
p lx p It 0 

To determine the first order terms of the flow properties behind the front, 
substitute of the expansion series (A) into equations (2.1) to (2.3) and use 

•jig 

the relation Up = Up. This gives three equations identical to equations 
(2.27) to (2.29) except that Up is changed to Up . Equations (5.17), 
(5.18), (5.19), and (5.22) together with these three equations form a system 
which can be shown to be redundant in a manner analogous to the behavior of 
equations (3.48) to (3.51). Thus, equation (5.21), for example, can be 
omitted. The remaining six equations form a linear system from which the six 
spatial and time derivatives of the flow quantities behind the reaction front 
can be calculated. The initial acceleration of the front is expressed 
directly by equation (5.19); it is independent of the piston acceleration. 


5.3 Solution of Flow Problems Involving a Detonation Process 

The detonation is assumed to be initiated by moving a piston into a 
nonuniform region (x > 0) with initial speed and acceleration u p and ap. 
Depending on the initial velocity of the piston, two types of detonation can 
occur. First, if u p is greater then or equal to the Chapman- Jouguet value 
Up discussed in the previous paragraph, the detonation is termed a strong 
detonation. (Fig. 7) 
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The flow relative to the reaction front is supersonic ahead and subsonic 
behind, which is the same as for an adiabatic shock. The differential 
equations of the flow field except at the reaction front are equations (2.1) 
to (2.3), and the jump conditions are equations (5.1) to (5.3). The only 
equation that is different from those describing the shock problem is the 
energy jump condition (5.3), which has an extra term to account for chemical 
energy due to burning. 

The procedure used in solving the problem is exactly the same as that 
for the shock problem. First rewrite equation (5.3) using equations (5.1) and 
(5.2) 

(Y ~ !)v 2 2 - (y + 1 ) v 1 v 2 + lc 2 + — (Y “ l)v 2 => 0 (5.23) 

V 2 “ v ! 

Now, substitution of the expansion series (A) into equation (5.23) and 
comparison of like powers of t gives, at the leading term 
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m 



which defines the initial speed of the reaction front. The first order term 
in t yields 



where E is defined by equation (2.33) with u g changed to u^. With equation 
(2.21) replaced by equation (5.25), and by making use of equation (5.24) for 
the initial velocity of the reaction front, the problem can be solved using 

exactly the same procedure described in Chapter 2. 

In the second case, in which Up is less than or equal to Up , a 
flow involving a Chapman-Jouguet detonation is possible. The flow immediately 
behind the Chapman-Jouguet front need not be equal to the piston velocity 
Up. Adjustment in this case is effected by a centered expansion wave which 
follows the detonation front immediately , since both the front and the first 
wave in the expansion fan move with sound velocity relative to the gas behind 

the front [9]. This expansion wave drops out only if Up * Up . The 
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problem is solved by combining the Chapman- Jouguet detonation and expansion 
wave results. First define the region ahead the detonation front as 
region II, the region behind the expansion fan as region I, and the fan itself 
as region III (Fig. 8). 



The procedure used is to solve the Chapman-Jouguet detonation with 
initial data in x > 0. Then the flow properties immediately behind the front 
are the initial conditions for the expansion fan problem. The solution in 
region II is completely specified by the initial data in x > 0. The 
Chapman-Jouguet detonation is solved with the results outlined in 
section 5.2. The leading terms of the flow properties immediately following 
the Chapman-Jouguet detonation will be denoted, p^^\ and They 

are all defined in equations (5.13) , (5.14), and (5.20) by replacing 
subscript ”1” by "4". The initial velocity of the front, u^, is defined by 
equation (5.16). Since these solutions are the initial conditions for the 
expansion fan problem, the leading terms of the flow properties inside the fan 
and behind the fan can be worked out using the equations in Chapter 3. 
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To determine the first order terms of the flow quantities behind the 

front, we note that an< * are now known f rom equations 

(5.17), (5.18), and (5.21). Because the front coincides with the leading 

characteristic of the expansion fan, of equation (5.12) is equal to 

x (£_,t) of equation (3.7a). Therefore 
c l 


and 



(5.26) 


(5.27) 


a^ is defined in equation (5.19), so k^j) is known. Now, rewriting the 
boundary condition (3.10) using (5.26), yields 


q 3 (l >(5 2 > 


'0 


(5.28) 


The quantities Q^u^ + Q^cq are those known from equations (5.17), (5.18), 
and (5.21). Using the expression (5.27) for and taking (5.28) as the 

boundary condition on OA, the first order terms of all the flow properties 
inside the fan and behind the fan can be calculated in the manner discussed in 
Chapter 3. 


5.4 Solution of Flow Problems Involving a Deflagration Process 

In deflagration processes the situation is in many respects quite 
different from that encountered with detonation processes. Suppose a weak 
deflagration wave begins at a piston at x = 0, t = 0. Then the velocity u 
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of the burned gas behind the deflagration front is negative, u < 0 [9], This 

is compatible with the conditions of the problem only if the piston is 
withdrawn with a speed at least equal to u. What actually happens is that a 
precompression shock is sent out into the explosive. It pushes the explosive 
gas ahead with a velocity just sufficient to ensure that it may attain the 
same velocity as the piston when it is swept over and burned by the 
deflagration front. The occurrence of a precompression shock is in complete 
agreement with, or rather a consequence of, Jouguet's rule that the flow ahead 
of a deflagration is subsonic. Consequently, the deflagration influences the 
state of the gas ahead of it. Define the region ahead of the precompression 
shock as region II, the region behind the shock as region V, and region 
between the piston and the Chapman-Jouguet deflagration as region I (Fig. 9). 


Particle c 0 t 
path 


Particle 

path 

i 



Figure 9.- Weak deflagration with precompression shock. 


To determine the leading terms of the flow properties in region I and in 
region V and the initial speed of the deflagration (u^) and of the shock (u 8 ), 
the results of the Chapman-Jouguet wave problem and the shock problem are 
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used. First, an equation which relates the initial velocity of the piston 


(u ) to the flow properties behind the shock is worked out. From 

r 

equation (2.24), two relations can be established for -u g ) an( ^ 

(u 2 ^^ - u g ) in terms of (u 2 ^^ - 



Using c = Yp/p from equation (2.25) and (2.26), 
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From equations (5.13) and (5.14), 
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Substituting equation (5.16) into (5.20) with the relation Uj^) » 
results in 


(5.30) 


(5.31) 


(5.32) 


(5.33) 
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(5.34) 


u„ i c,(«> - U J°> » c,-0) 


1 ♦ ail.. : »> 


L\l 


2c, 


( 0 )‘ 




q ( y - i) 

* 6 “' 2 


Equation (5.34) expresses the desired relation. The solution for u^®^ can 
be determined iteratively by using equations (5.30) to (5.34). Once u^^ 
is determined, the velocity of shock, the leading terms of the flow properties 
behind the shock, the velocity of deflagration and the leading terms of the 
flow properties behind the deflagration can be calculated successively using 
the appropriate results from earlier chapters. 

The first order terms of the flow properties in region V and region I, 
and the initial accelerations of the Chapman-Jouguet reaction front (a^) and 
of compression shock (a Q ), consist of 14 unknowns to be determined. The 
equations to be used are (i) the three jump conditions across the compression 
shock, equations (2.19), (2.20) (with region I there changed to region V), and 
equation (5.25); (ii) the conditions across the Chapman-Jouguet reaction 
front, equations (5.17), (5.18), and (5.19) (changing subscript 2 to 5 and and 
using the lower sign for deflagration); (iii) equations (2.27) to (2.29) for 
conditions inside region I; (iv) the counterparts of equations (2.27) to 
(2.29) for region 5, which are 


( 0 ) . ( 0 ) 
p 5x u 5 + p 5t c 0 ' -p 5 u 5x 
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(v) the boundary condition a p * u i x u p + u lt c 0» an< * equation (5.21) which 

follows from the fact that Mach number is unity behind the Chapman- Jouguet 
reaction front. The 14 unknowns can be determined by solving these 14 
algebraic equations. 

It is of interest to note that since the Mach number behind the Chapman- 
Jouguet reaction front is unity, this front is a characteristic line. 

It was concluded in Chapter 3 that the initial acceleration of the piston, a p , 
does not affect the flow properties ahead of the tail characteristic of an 
expansion fan, and does not affect the acceleration of the tail characteristic 
line itself. Since the Chapman-Jouguet reaction front is a characteristic 
line, the same situation would be expected to be true here; that is, the flow 
properties in region V and a g and a^ should be determinable without 
knowing a p . It is not obvious from the above set of 14 equations that this 
is so. It can be seen to be true by the following argument. Consider the 14 
algebraic equations listed above, omit ill (equations 2.27 to 2.29) and 
v (a p ■ u ix u p + u lt c o^* an< * i nc lud e t * le Riemann invariant relation (3.46) for 
the Chapman-Jouguet deflagration front. This results in 11 equations, which 
can be used to solve for the 11 unknowns, p^, P 5 t > Pg x > P 5 t » U 5 X * u 5t» a s* 
a^, Pi^> Pi^> an< * u i^* these, all quantities in region V and 

are known. Thus, the present case is consistent with previous result that 
the piston acceleration only affects the flow properties behind the tail 
characteristic, but not the initial acceleration of the tail characteristic, 
or the flow properties ahead of the tail characteristic. 

This completes the discussion of the elementary piston problems with 
chemical reaction. First, the Chapman-Jouguet solution has been worked out in 
detail in section 1. The strong detonation problem, which is similar to the 
adiabatic shock problem, is discussed next. The Chapman-Jouguet detonation 
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with an expansion fan following the detonation front and the Chapman-Jouguet 
deflagration with a precompression shock are discussed in sections 3 and 4, 
respectively. In these problems, explicit solutions are not worked out in 
detail because each involves extensive algebra; however, they can be 
determined numerically easily using the listed linear systems of equations. 
Once the separate elementary problems have been solved, they can be combined 
with the shock or expansion fan solutions in the preceeding chapters using 
exactly the same procedure as in chapter 4 to solve initial value problems 
with chemical reaction. Two examples of such problems are indicated in 
figures 10 and 11. Figure 10 shows a typical wave pattern including chemical 
reaction and an adiabatic expansion fan; Figure 11 illustrates a pattern with 

chemical reaction and a shock wave. 

Pressure at t ■ tj 



Particle t line 



Figure 10.- Wave pattern including chemical reaction and 
adiabatic expansion fan. 
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Figure 11.- Wave pattern with chemical reaction and shock wave. 


V 


* 


65 



CHAPTER 6 


CONCLUDING REMARKS 


The purpose of this work is to discuss the initial value problem of one- 
dimensional gas-dynamics involving discontinuous, nonuniform initial data. 
Canonical solutions which are valid in a small x, t region around a dis- 
continuity, and which include the first order effects of nonuniformities, have 
been derived explicitly. Similar solutions corresponding to discontinuous but 
otherwise uniform initial data have been used in the past as building blocks 
for the numerical solution of general initial value problems. The results 
presented in the present work, in which the nonuniformities in the data are 
explicitly taken into account, are intended to be used as building blocks in 
an improved numerical scheme which will permit the use of much larger mesh 
sizes than do previous methods. 

The theory has been derived by considering a group of elementary piston 
problems. Solutions with a shock or with a centered expansion wave have been 
worked out individually in order to relate initial flow properties and their 
gradients to the speed and acceleration of the discontinuity waves. Then they 
have been combined to represent the solution of an initial value problem by 
regarding the piston path as a contact line. In addition, problems with 
chemical reaction are discussed in terms of elementary piston problems which 



involve strong detonation waves, Chapman-Jouguet detonation waves, and 
deflagration waves. These solutions can be combined with the shock or 
expansion fan results to derive canonical solutions of the initial value 
problem with chemical reaction. 

The actual implementation of the current results into an improved 
numerical scheme for general initial data is beyond the scope of the current 


paper. 
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Jouguet detonation waves, and deflagration waves. 


17. Key Words (Suggested by Author(s)) 

Unsteady flow 


18. Distribution Statement 

Unclassified - Unlimited 

Initial value problem 
shock wave 
Chemical reaction 


Subject Category 34 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

UNCLASSIFIED 

UNCLASSIFIED 


73 

A04 


N-305 


For sale by the National Technical Information Service, Springfield, Virginia 22161 



















